Mon. Not. R. Astron. Soc. 000, 000-000 (2010) Printed 28 January 201 1 (MN ETeX style file v2.2) 
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ABSTRACT 

The peculiar motions of galaxies can be used to infer the distribution of matter in the 
Universe. It has recently been shown that measurements of the peculiar velocity field indicates 
an anomalously high bulk flow of galaxies in our local volume. In this paper we find the 
implications of the high bulk flow for the power spectrum of density fluctuations. We find 
that analyzing only the dipole moment of the velocity field yields an average power spectrum 
amplitude which is indeed much higher than the ACDM value. However, by also including 
shear and octupole moments of the velocity field, and marginalizing over possible values for 
the growth rate, an average power spectrum amplitude which is consistent with ACDM is 
recovered. We attempt to infer the shape of the matter power spectrum from moments of the 
velocity field, and find a slight excess of power on scales ~ 1ft, -1 Gpc. 

Key words: cosmology: large scale structure of the universe - cosmology: observation - 
cosmology: theory - galaxies: kinematics and dynamics - galaxies: statistics 



; 1 INTRODUCTION 



Peculiar velocities are useful cosmological probes. In principle, the 
peculiar velocity field is an unbiased tracer of the underlying matter 
distribution, and should be sensitive to structure on scales larger 
than the nominal size of the survey. The study of large scale flows 
has a long and hallowed tradition; throughout the 1980s and early 
1990s, measurements of large scale flows were deployed to not only 
constrain the fractional energy density in matter, fl m but also to 
identify possible sources of gra vitational attraction which might lie 
outside current galaxy surveys tetrauss&WiiiIckUl995l) ). 

Velocity catalogues are hard to construct. While current red- 
shift surveys have of order 10 5 objects, velocity catalogues are 
typically restricted to of order 10 galaxies. Yet over the past 
three decades, enough measurements of peculiar velocities have 
been accrued to be able to construct a reasonably complete cat- 
alogue out to a maxim um distance of about 100/i~ Mpc. In a 
recent series of papers dWatkins et al"1l2009l : iFeldman et aUl2010h . 
the authors showed that a collection of peculiar velocity surveys 
could be combined to construct a reliable, well-behaved "compos- 
ite" catalogue. This composite catalogue was then used to extract 
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the simplest statistic: the bulk flow. The results were surprising: 
there is clear evidence for a large bulk flow, which is very unlikely 
within the current preferred model of structure formation, i.e. a flat, 
Friedman-Robertson- Walker universe with an appreciable cosmo- 
logical constant permeated by Gaussian, scale-inva riant perturba- 
tions ( k nown as the ACDM m odel). The findings of Watki ns et al.l 
(2009); Fel dman et all feoicli leads to the question: is there excess 
power on large scales? As yet we do not have a direct measurement 
of the power spectrum of density fluctuations on these scales in the 
nearby Universe; the cosmic microwave background does probe a 
wide range of scales but at high redshift. Could the large bulk flow 
be a signature for large scale fluctuations at low redshift? It is this 
question we wish to address by using the measurement of the bulk 
flow and a few of the other lower moments to estimate the power 
spectrum of density fluctuations. 



|jaffe& Kaiser! d 1995b inferred constraints on the matter power 
spectr um parameterized by as and shape parameter T with data 
from [L auer & P ostman! a 19941) . and seperately with data from 
iRiess e t al. ( 1995). They found that including only the bulk flow of 
the Lauer & Postman catalogue favoured a model with high large 
scale power, which became more reasonable when the shear of the 
velocity field was included, suggesting a very large scale excess of 
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power. However, the data from lRiess et alj i 19951) did not suggest 
any large scale excess. 

The shear and dipole moments used bv ljaffe&Kaisen l ll995h 
are slightly different to the minim um variance shear an d dipole 
moments we use here, because the Ijaffe & Kaisei] d 19951) method 
assumes that the non-modeled ve l ocity is noise. Thus the dipole 
moment from the Jaffe & Kaisei] J 19951) method will depend on 
whether the shear is included, whereas the minimum variance for- 
malism we employ here estimates each moment inde pendently of 
the to tal number of moments analyzed. Further, Ijaffe & Raised 
( 1995) compare very different catalogues with different geome- 
tries, i.e., densities, distribution, etc. The minimum variance 
method combines individual measurements to estimate the velocity 
moments for an idealized geometry, and the results are thus inde- 
pend ent of the geom e try of any particular survey. 

iKolatt & Dekel J 19951) inferred constraints on large scale 
structure from the POTENT velocity field. Uncertainty in the 
growth of structure rate, /, was included by reporting the degen- 
erate combination of P(k)f 2 . With current improved constraints 
on the growth rate, and more peculiar velocity data available, we 
attempt to constrain the unbiased matter power spectrum from mo- 
ments of the peculiar velocity fi eld. 

Recentlv lSong et all j201Ch applied a method to estimate bulk 
flows from two point correlations in galaxies distributions, and 
found no evidence for a large bulk flow in data from SDSS. How- 
ever, their analysis only probes scales up to 0.03 /iMpc -1 , whereas 
we see later that the anomalously high bulk flow is most sensitive 
to scales >0.01 /iMpc -1 . 

This paper is structured as follows. In section[2]we review the 
methodology to incorporate an individual peculiar velocity mea- 
surement into an esti mate of our loc al peculiar velocity field. We 
review the method of iKaiseJ Jl988h to weight peculiar velocity 
m easurements to estima te the bulk flow. We also review the method 
of IWatkins et al. I J2009I) to use peculiar velocity measurements to 
estimate the bulk flow of our local volume, so that results from dif- 
ferent surveys can be directly co mpared. We describe the peculiar 
velocity catalogu es compiled bylWatkins et al] (2009). We then re- 
view the work of Feldman et al. (2010) to include higher moments 
of the velocity field. In section [3] we describe how we relate these 
moments of the velocity field to constraints on large scale structure. 
We present our results in section|4] 



2 COMBINING PECULIAR VELOCITY 
MEASUREMENTS 

To estimate the peculiar velocity of galaxies at cosmological dis- 
tances, we need both the measured redshift of the galaxy, and an in- 
dependent distance measure, such as the luminosity distance. From 
the distance measure, we can estimate what we would expect the 
redshift to be solely due to the Hubble flow. We can then attribute 
the difference between this expected redshift and the measu red red- 
shift to the peculiar velocity of the galaxy JPeeblej J1993I) ). There 
are two key difficulties with the method. The first is that the uncer- 
tainty on the luminosity distance is typically rather large: ~ 10 to 
20% for Tully-Fisher, Faber- Jackson, or Fundamental Plane mea- 
surements, and ~ 5% for supernovae. The second difficulty is more 
inherent: this method only provides the line of sight component of 
the peculiar velocity. The importance of this effect can be illustrated 
if we consider a hypothetical bulk flow from an arbitrary north to 
the south. A survey of galaxies to our north or south would be sen- 
sitive to this flow, whereas a survey to our east or west would not. 



The approach we describe here is to combine individual peculiar 
velocity measurements to estimate the velocity field of our local 
volume by weighing each measurement according to how sensitive 
it is to each component of the underlying velocity field. We can de- 
scribe the velocity field v(r), where v is the three components of 
the peculiar velocity field at position r. For each galaxy (labelled 
n), we measure the line of sight component of this field, S, i.e. 
S„(r) = v(r n ) • f„ where f n is the unit vector pointing in the 
r n direction. We assume this measurement is drawn from a Gaus- 
sian distribution with variance a n . An additional term, a 2 , is in- 
cluded in the variance, to account for non linear flows. The method 
is fairly insensitive to the particular choice of <r„ , as the combined 
uncertainty of a^ + a 2 tends to be dominated by the measurement 
uncertainty, a\. 

The simplest result we can quote for a survey of peculiar ve- 
locities is the 'bulk flow' vector, u, the average velocity magnitude 
and direction of the galaxies in a survey. We must be careful when 
combining peculiar velocity measurements to include the effect of 
only measuring the line of sight component. This is achieved by 
multiplying each component of S n (in the coordinate system of the 
bulk flow) by a weight w a , which essentially depends on the or- 
thogonality of S to the coordinate system of the bulk flow, so that 
u a = Yin w " "*5"- lK a i se d < fl988l) showed that the weights for the 
bulk flow are 
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Because of sparse sampling, bulk flow results between 
different surveys ar e not necessarily immediately comparable. 
Watki ns et al. U2009I) devised a method to weight the bulk flow of a 
particular survey to estimate the bulk flow of a hypothetical spher- 
ically symmetric survey, where the density of galaxies n falls of 
as n(r) oc exp(— r 2 /2R 2 ), where R 2 is the characteristic depth 
of the survey. This allows bulk flows from different surveys to be 
directly compared, and combined to give us a better estimate of 
the bulk flow of our l ocal v olume of space. We start with the for- 
malism from [Raised dl988l) . and then introduce the method from 
IWatkins et al.U2009l) . 

The idea is that the bulk flow we measure, u, is essentially a 
convolution of the peculiar velocity field v with the window func- 
tion of the survey W, 



Ui(r ) — J d 3 rWij(r)vj(r + r ) 



(3) 



This allows us to calculate the variance of u directly from the con- 
volution theorem 



(muj) = / d 3 kP v (k)W 2 (k) (4) 
where the window function is given (in real space) by 
Wa(r) = A^ , £S{r-r n )J^- s (5) 



Rjj , due to measurement variance, and a noise term, R) 1 - 1 due to 

the variance from non linear flows, so that (iliUj) = R^ + R t j ■ 
We then have that 

R (e) - A- 1 



It is extremely useful to split {uiUj) into two terms: one term, 



?C<0 



(6) 



and 

= J d 3 kW?j(k)P v (k) (7) 
where 

W£ = Wa{k)W; m (k)hkm (8) 

We will return to equations[6]and[7jwhen we come to estimate large 

scale structure in section|3] 

The approach of Watki ns et al, I d2009l) is to find weights that 
minimize the average variance between a particular bulk flow vec- 
tor u and the average bulk flow of our hypothetical survey, U, 
we assume the measurement error is uncorrelated 
with the bulk flow, we can expand ((u^ — Ui) 2 ) as 

((Ui - Ui) 2 ) = w i,n w i,™ (SnS m ) + (U 2 )-2^2w itn (UiV n ) 

Before minimizing this expression with respect to Wi t „, we enforce 
the constraint that W 2 j (k) — > 1/3 as k — > 0, so that the weighting 
is equal for each dimension of the peculiar velocity field. We thus 
have to minimize 

^ Wi : „W i:m {S n Sm) + {Ui) 
n,m 

-~2^2,w lt n (uiUi) \(P„ m Wi,nWi,m - t:) (10) 



where A is a Lagrange multiplier and 




Differentiating equationQI)] and setting the result equal to zero 

gives 



I'll )Wi,m = (S n Ui) (12) 

m 

It is easier to solve for the weights, tOj, m , if we rewrite equation 
I12lin matrix form. Substituting G for {S„S m ), P for P nm , Wj for 
iOi, m , and Q ; for (S n Ui) gives 

(G + AP)w i = Qi (13) 

which is now easy to solve for the weights Wi. 

wi = (G + AP^Qi (14) 

As well as the bulk flow of the survey, we can also consider 
higher moments of the peculiar velocity field if we consider the 
field as a Taylor expansion, as in 

Vi(v) = Ui + Uijrj + Uijkrjr k + . . . (15) 

Ui is the bulk flow vector - also known as the 'dipole moment'. 
Uij is the 'shear tensor' or the 'quadrupole moment', and provides 
information about the distance at which the bulk flow attractor 
is located. The 'octupole tensor' Uijk provides inf ormation about 
the ve locity field on scales smalle r than the survey. iFeldman et ail 
(2010) have extended the work of W atkins et alj d2009j) to include 
these higher order moments. If we assume that the peculiar veloc- 
ity field is entirely due to gravitational infall, we expect the field to 
be curl-free, and consequently Uij and Uijk to be symmetric. Thus 
the peculiar velocity field can be described to third order by the 19 
independent velocity moments, m . 

Before we proceed, we have to be careful with the definition of 
the octupole tensor, because some components of the tensor overlap 
with the dipole moment. As such, we modify the expansion of the 
velocity field to be 
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Wi(r) = Ui + Uijrj + Uijk {r 3 r k - A jfe ) + ■ ■ ■ (16) 
where Ajk is given by 

A 3fc = / r jrk d\ (17) 
Jv 

in order to remove overlapping components. The line of sight com- 
ponent is then 

s(r) = Uifi + Uijrfifj + Uijk (r 2 fifjf k - A jfc fi) + ■■■ (18) 

Feld man et al] l l2010h have calculated the minimum variance 
weights for the 19 components of the third order expansion, over k 
ranges from 0.002 to 0.196 /iMpc~\ We can think of the resulting 
set of 19 components of the dipole, quadropole and octupole as a 
form of data compression containing the highest signal to noise in- 
formation in a given peculiar velocity survey. Indeed, that will be 
our philosophy - to use this form of the data to infer information 
about the underlying density field. 

We stu dy the 'COMPOSIT E' peculiar velo city catalogue 
compi led bv lWatkins et al.l ( l2009t) . and also used in lFeldman et al.l 
(2010). The catalogue consists of 4,536 peculiar velocity measure- 
ments, with a characteristic depth of 34 /i _1 Mpc. The characteristic 
depth is given by the average distance to each galaxy, weighted by 
the inverse square of the peculiar velocity uncertainty. 



3 RELATING VELOCITY TO MATTER 

We have presented a method to combine line of sight estimates of 
the peculiar velocity of individual galaxies into an estimate of the 
moments describing the velocity field of our local volume. We now 
want to be able to compare this measurement to expectations from 
our cosmology. The basic idea is that, in the linear regime, galaxies 
flow towards local over-densities of matter, so the velocity field 
v(r) is given by 

v(r) = ^/dV5(r')^4 (19) 
47T J | r ' - r| 

/ is the perturbation growth rate, d In S/d In a, and S is the matter 
density contrast. The matter density contrast is modelled as a Gaus- 
sian random field with a power spectrum defined as (|<5(fc)| 2 ) = 
(27r) 3 P(fc), where 8 is the Fourier transform of the density con- 
trast in real space. 

A useful way of parametrizing P(k) is in terms of band pow- 
ers: 

P ^ = { P otherwise ^ ™ 

To estimate the most likely matter power spectrum, based on 
the peculiar velocity data, our approach is to construct, and then 
minimize, the likelihood function £ as a function of P a : 

-2\n£(P a ) ocln|C| + u T C _1 u (21) 

where u are the velocity moment components. The covariance ma- 
trix, C, is derived from equation|4] 

Cpq = Rpq + Ppq (22) 

where the 'error matrix' R p E q is given by equation [6] and the 've- 
locity matrix' R p V q is given by equation [7] We relate the velocity 
covariance matrix to the power spectrum by equation [7] so that 
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Figure 1. Likelihood distributions for the average P(k) value in the k range 
of our window function, inferred from a simulated peculiar velocity cata- 
logue. Whether we include only the bulk flow of the velocity field (labelled 
'Dipole'), or also include the shear and octupole moments (labelled 'Oc- 
tupole'), we recover a value which is compatible with the fiducial ACDM 
value. For clarity, the likelihood functions have been scaled arbitrarily to 
the maximum value of the dipole function. The ACDM value has been il- 
lustrated here and in further such plots as a Gaussian distribution with an 
arbitrary standard deviation of 10% of the fiducial value. 

We can choose the width of each band-power by integrating 
the window function (which is independent of P a ) over the bin 
range k a to fc Q +i. We can then factor in the average power ampli- 
tude in this range, so that 

f 2 

a 

where the kernel /C is given by 

= [ ^ dkW* q (k) (25) 

J k a 

We thus have a likelihood function for the velocity moments in 
terms of the growth rate and power spectrum band-powers. 

We have tested our methodology on simulated peculiar veloc- 
ity catalogues, generated from a ACDM power spectrum. We gen- 
erate a series of mock catalogues from a set of N-body simulations, 
estimate the corresponding dipole, quadropoles and octopoles and 
then minimize the corresponding likelihood to recover the input 
power spectrum. We illustrate our results for one such realization 
in Figures Q] and [2] and confirm that our method does not generate 
any spurious large scale power. 



4 RESULTS AND DISCUSSION 

We first consider the most likely total average power of our sam- 
ple, over the entire k range of our window function, k — 0.002 to 
0.196/i Mpc" 1 . We see in Figure [3(a)] that when we include only 
the dipole moment of the velocity field, the average power for each 
survey is much higher than the A CDM value. This is e ntirely con- 
sistent with the results found in Wat kins et al. I d2009l) . When we 
also include the shear and octupole moments, we find much better 
agreement with the ACDM value, as we can see in Figure [3(b)l This 
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Figure 2. In the upper figure is plotted the power spectrum shape, inferred 
only from peculiar velocity data generated from a simulated catalogue. We 
have chosen three bands, evenly spaced in log k. The shaded region rep- 
resents the marginalized lcr uncertainty on the power in each band. For 
comparison, the theoretical ACDM power spectrum has been overplotted. 
The results shown here include the dipole, shear, and octupole moments 
of the simulated velocity field. As discussed later, the uncertainties in each 
band have also been marginalized over the growth rate. In the lower panel 
are plotted the average values of the dipole, shear, and octupole window 
functions. The dipole is most sensitive on scales k <0.01 hMpc" 1 . 



is similar to the effect observed bv |jaffe&KaiseiUl995h when in- 
cluding the shear of t he peculiar velocity field, and also noticed by 
lFeldmanetaT] ( l201(]h . 

In Figure [3] the growth rate was fixed at the fiducial ACDM 
value, to illustrate the importance of including higher moments of 
the velocity field. However, presentl y the best con s traints on the 
growth rate at low redshift are from iPeacock et al. I fcOOlIV s mea- 
surement of the redshift space distortion compression parameter. 
This constrains the growth rate to / = 0.49 ± 0.14. We now con- 
sider the effect of marginalizing over the growth rate, with this 
prior applied. Likelihood contours for the growth rate and average 
power in the COMPOSITE survey are shown in Figure|4] When we 
marginalize over the growth rate, we see in Figure [5] that the aver- 
age power does not exclude the fiducial ACDM value. However in 
the context of ACDM models, the growth rate is not a free parame- 
ter, rather it is determined primarily by Q. m (which also affects the 
shape of the power spectrum and the average power). In order to 
obtain consistency with the predicted average power, on requires 
/ ~ 0.7, which is much larger than the fiducial ACDM value of 
0.48. 

In addition to the total average value, we can attempt to infer 
the shape of the power spectrum. The window function was divided 
into three ranges, evenly spaced in log k. The likelihood function 
was then evaluated in terms of the power in each of the three bins, 
P a , and the growth rate. As before, a prior on the growth rate was 
applied, and this parameter was marginalized. When we only con- 
sider the bulk flow, we find that the likelihood is sufficiently broad 
and correlated that it is, in practice, impossible to pin down three 
independent estimates of P a . This is disappointing - the measure- 
ment of the dipole alone does not allow us to identify if there is 
large scale power. 
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Figure 3. In the upper panel is the likelihood distribution for the average 
amplitude of the power spectrum, inferred only from the dipole moment. 
To account for the anomalously high bulk flow, we have to conclude a 
power spectrum amplitude which is incompatibly higher than the ACDM 
value. In the lower panel, we show the effect of including the shear and oc- 
tupole moments of the velocity field. By now including the higher moments, 
the results are generally more compatible with the ACDM value. Here, the 
growth rate, /, is kept fixed at the fiducial ACDM value. In both figures, the 
shaded area is the lcr confidence region. The shaded area of the upper figure 
appears disproportionately wide because the large tail of the distribution is 
distorted by the logarithmic scale on the horizontal axis. 

However, when we also include the shear and octupole mo- 
ments of the velocity field, we find that the band powers are vir- 
tually uncorrelated, and we can obtain three independent measure- 
ments. In Figure[6]we plot the marginalized la uncertainty for each 
P a for the COMPOSITE catalogue. There is a slight detection of 
excess power on the largest scales. 

We find that our estimate of the excess power seems to re- 
flect the large scale estimates of P(k) found in the estimate from 
both t he SDSS main galaxy power spectrum from Percival et al. 
(2007) and the 2dFGRS power spectrum from lCole et al.l fc005). 
One should be careful about over-interpreting these similarities, as 
both the SDSS and 2dFGRS estimates of power on large scales 
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Figure 4. Likelihood contours for the COMPOSITE survey, with dipole, 
shear, and octupole moments included, with contours spaced at 1, 2 and 
3cr confidence levels. We are now considering the average power and the 
perturbation growth rate as free parameters. A prior on the growth rate has 
been applied. For comparison, the 1 dimensional likelihood distribution for 
the COMPOSITE survey in Figure [3(b)] is a slice through this likelihood 
contour at the fiducial ACDM value of the growth rate at z = 0, i.e. / = 
0.48. As can be seen here, including observationally constrained values of 
the growth rate allows for a wider range of allowed P a . 
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Figure 5. Average power for the COMPOSITE survey, with dipole, shear, 
and octupole moments included, now marginalized over the growth rate. 
The shaded region indicates the lcr uncertainty boundaries on the likelihood 
function. The expected ACDM value is now included in this range. 

are, in principle, severely affecte d by edge effects. Indeed , the re- 
cent analysis of the SDSS data in lSvlos Labini et al] ^2009) should 
make one be wary about over interpreting estimates of structure on 
the largest scales of the survey. Yet the fact that the same level of 
fluctuations is obtained from a peculiar velocity survey, subjected 
to a completely different analysis may be an indication that redshift 
survey estimates of structure on the largest scale must be taken se- 
riously. Furthermore, an excess of cl ustering on similar sc ales at 
higher redshift was recently found by iThomas et alj d2010|) in the 
MegaZ DR7 photometric redshift survey. 
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Figure 6. Power spectrum shape inferred from dipole, shear, and octupole 
moments for the COMPOSITE catalogue. The shaded regions are la un- 
certainties, marginalized over the growth rate and the other P a . We find a 
slight excess of power on scales of ~ Ih^ 1 Gpc. These are the scales at 
which we are most sensitive. This excess of power appears to agree with 
the lar gest scales in the 2d FGRS and SDSS surveys, which are plotted here 
as in IPercival et all llOOH) and have been deconvolved from their survey 
window functions. 



What should we make of this excess of power, and how does 
it relate to our initial question: the anomalously high bulk flow? To 
start with, we have demonstrated that it is incomplete to infer re- 
sults from measurements of the bulk flow alone, as including higher 
moments of the velocity field yie lds significantly different results. 
As noted in iFeldman et al] (2010), the low shear moments of the 
velocity field suggest that the sources responsible for the bulk flow 
are very far away. In the formalism we have adopted here, only 
one of the bulk flow moments is anomalously high - the other bulk 
flow moments, and the shear and octupole moments, are consistent 
with ACDM. From one perspective one could argue that the extra 
freedom from including 19 moments of the velocity field, and the 
growth rate, as opposed to just three moments of the dipole mo- 
ment, provide a way to interpret the anomalously high bulk flow 
which is consistent with ACDM. On the other hand, all moments 
are not equivalent. The dipole moments are special in the sense 
that they probe the largest scales (and are also the moments that are 
most robustly measured). Consequently the inferred shape of the 
power spectrum is different from ACDM. 
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